Casimir effect with dynamical matter on thin mirrors 
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We calculate the Casimir energy for scalar and gauge fields in interaction with zero-width mirrors, 
including quantum effects due to the matter fields inside the mirrors. We consider models where 
those fields are either scalar or fermionic, obtaining general expressions for the energy as a function 
iq , of the vacuum field fPI function. We also study, within the frame of a concrete model, the role of 

the dissipation induced by those degrees of freedom, showing that, after integration of the matter 
fields, the effective theory for the electromagnetic field contains modes with complex energies. As 
for the case of Lifshitz formula, we show that the formal result obtained by neglecting dissipation 
coincides with the correct result that comes from the quantum fluctuations of both bulk and matter 
Q_l' fields. 
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Cn ' It has become increasingly clear that the use of sensible models for the description of the media that 

constitute the mirrors is an unavoidable step for the refinement of Casimir energy calculations [l|. In 
particular, the theory of quantum open systems is the natural approach to clarify the role of dissipation 
in Casimir physics. In this framework, dissipation and noise appear in the effective theory of the relevant 

(T^ degrees of freedom (the electromagnetic field) after integration of the matter degrees of freedom. 

In this letter we consider mirrors described by thin films, equipped with matter fields confined to them; 
these localized fields shall induce an effective action for the coupling of the vacuum fields to the mirrors. It 

qq ' will turn out that the resulting effective action falls into a class of model that we had studied in a previous 

f^ , work [2[ . We then apply those results here in a quite straightforward way, obtaining the Casimir energy for 

different cases: regarding the vacuum fields, we shall consider the cases of a massless (real) scalar field and 
of an Abelian gauge field. As we shall see, under our general assumptions, the Casimir energy of the latter 
may be obtained as a derived result of the former, after performing some identifications. 

The Casimir effect produced by thin plasma sheets has been previously considered by other authors 
(see [3, |j] and references therein) imposing appropriate boundary conditions on the position of the sheets [3[ , 
considering a fluid of non relativistic electrons confined to the surface [4|], or describing the matter inside 
thin films using the particle in a box model Q . We will extend these results to the case of relativistic bosonic 
or fermionic degrees of freedom confined to the mirrors. Besides, for a particular fermionic model, we will 
show that, in the effective theory for the electromagnetic field, the modes have complex energies, due to the 
possibility of transferring energy from the bulk to the thin surfaces. The Casimir energy is of course always 
real, and given by the sum of cigcnfrequcncies of the full system (bulk and matter fields). 

Let us first consider the effective action, S^Jf, for a real scalar field <p, which results from its interaction 
with matter confined to two identical, flat, parallel and zero-width mirrors in d spatial dimensions, defined 
by the equations Xd = and Xd = a. The total Euclidean action S for tp will then be of the form: 

s(ip) = s (<p) + s<3m , (1) 
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where So = \ J d d+1 x d^ipd^p, defines the free theory. S e Jf do, of course, depends on the nature of the fields 

confined to the mirrors and on the structure of their interaction with tp. The functional form of S e J. will be 
assumed to be quadratic in ip. This approximation is justified by the following reason: since the media are 



assumed to impose (approximate) Dirichlet boundary conditions, the relevant configurations correspond to 
small values for the fields on the mirrors (otherwise there would not be approximate boundary conditions 
there). It is then legitimate to assume that higher-order terms in those fields shall be strongly suppressed 
in comparison with the lowest, non trivial one (quadratic), since they will have a smaller weight in the 
functional integral. 

With this remark in mind, taking into account the fact that the matter fields are confined to the mirrors, 
and recalling that the interaction is assumed to be local, we have: 



S$A<p) = S$A<P) + S$A<P) (2) 



where: 



Sif f (<p) = \ I dx I dx' I d^sn / d d - x x\ 

x <p(xa,-x.\\,Xd)\{xQ - x' Q \-x.\\ -xj|)<5(.T d - a a )5(x' d - a a )ip(x' ,x.'px d ) , (3) 

with a\ = and a 2 = a. X|| denotes the d— 1 coordinates parallel to the mirror: x\, X2, ■ ■ ■ , x d —i- 

The interaction term will be obtained as an effective action coming from the integration of the degrees 

of freedom confined to the walls that interact locally with the scalar field ip(x) at x d = and x d = a. To 

simplify the notation in the following calculations, we shall denote by y a the spacetime coordinates for the 

matter fields, with y a = x a , for a = 0, 1, . . . , d — 1. 

Using a Fourier transformation with respect to the y a coordinates, 

5 S(V) = \j^^(k,a a )\(k)^k,a a ) (4) 

where k = (fc a ) a Z , are the Fourier space variables associated to the mirror's spacetime coordinates. Thus, 
the effect of the matter modes can, for this case, be encoded in a single function X(k). Moreover, depending 
on whether the dynamics of the matter fields is assumed to be relativistic or non relativistic, it will depend 
on k = y/k a k a , or on fco and |k|| |, respectively 

We see that the term ((4]) is the e — ► version of the models we have introduced in [4]. Thus, the result of 
the derivation presented there for the Casimir energy has direct applicability here. That is true, regardless 
of the microscopic model assumed for the matter field; however, before writing the expression for the energy, 
let us derive an expression for A from a concrete model. Note that, since the structure of S^ll depends 
trivially on the value of a a , we just need to consider the a = 1 case, shifting the field afterwards. 

A first model emerges from the assumption that S eff = S e ff is due to the linear coupling of ip to a 
microscopic real scalar field £(y), in d spacetime dimensions, equipped with an action S rn (^): 

where g is a coupling constant. The matter field £ may have a self-interaction, controlled by an independent 
coupling constant, implicit in S m . 

To proceed, we denote by W(J) the generating functional of connected correlation functions of £ (in a 
rf-dimensional spacetime), related to Z(J), the one for the full correlation functions: 

Z(J)= f De e -^(0 + /^(v)«v) _ (6 ) 
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by W = InZ. We then have that S e ff(ip) = — W[ig(p(y, 0)]. On the other hand, since only the quadratic 
part will be retained, 

S eff (<p) = -W[ig<p(y,0)] 

~ \d 2 J d d y J d d y' v(y,0)WW(y,y>My',0) , (7) 

where W^ 2 ' is the connected 2-point function. Finally, since we assume translation invariance along the 
mirrors' surface, we conclude that: 

A(k) = 9 2 W {2) {k). (8) 

Thus, for this first example, A may be read from the full propagator for the material field on the mirror. In 
the relativistic case, we know that: 

where z and II denote the wave function renormalization function and the full self-energy, respectively. Of 
course, we assume that the 2-point function has been renormalizcd, so that there is no remaining ambiguity 
in either one of these functions, which depend only on the d-scalar k. Since the full propagator may be 
related (via analytic continuation) to response functions, we see that the renormalization conditions (when 
needed) can be naturally identified with the static response functions for the field, due to the microscopic 
degrees of freedom on the mirrors. 

Thus, unless quantum effects introduced qualitatively important changes in the form of the full propagator, 
the conclusion is that, in this example, one should expect to have a good description by using the free 
propagator, namely, 

A M ~ -2^ > ( 10 ) 

with the obvious changes if the matter field were described by a different dispersion relation. 
The zero-point energy density £q can be written as 

£0 = Jta^lB(f), (11) 

where 

(12) 

This expression for the vacuum energy is of course equivalent to the sum over eigenfrequencies of the full 
quantum system. For the case of two mirrors, one is interested not in £ , but rather in the subtracted 
quantity, £q, defined as the difference 

£ (a) = £ o - £o(oo) (13) 

where £o(oo) denotes the surface energy density when the mirrors are separated by an infinite distance. 
Using the results of Ref . [2j , we can write the general expression for the Casimir energy: 
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The standard behaviour for A (re) is ~ n~ 2 for large values of re; this implies that the UV behaviour shall 
be softened with respect to the perfect mirror case, as expected. The IR region of the integrand is, on the 
other hand, strongly affected by the inclusion of the 2-point function; this is more evident for the case of a 
massless microscopic field. 

Let us now consider the case of an Abelian gauge field A^ in d + 1 dimensions, whose free action is of the 
usual Maxwell type 

So(A) = Jd d+l x^ u F^. (15) 

The coupling to the charged fields on the mirror, whose (total) electric current is J a (y) will be assumed to 
be of the standard J a A a type, where the current component Jd does not appear (we assume that no charge 
can escape from the mirror). Under these assumptions, the effective action can be written in the following 
fashion: 

e- s «fM = ±±^ = ^^ , (16) 

where we use "/j," to denote an unspecified field (or fields) which carry a conserved current, and e is the 
coupling constant. It is straightforward to verify that this defines a gauge- invariant functional of A a : 
S e ff(A) = S e ff(A + dui), for every smooth function lu, as a consequence of current conservation. 
Then, in the quadratic approximation for S e f/(A), 

S eff (A) ~ \j 2 J d d y J d d y> A a (y,0)li ab (y,y') A b (y',0) , (17) 

the kernel H ab is given by the current-current correlation function: 

n ab (y,y') = (J a (y,o)j b ( y ',o)) , (18) 

where the average symbols correspond to matter fields functional averaging, 

fVn...e- s ™(^ 
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Gauge invariance implies that H ab is transverse; this is more easily formulated in Fourier space: k a Tl ab (k) = 0. 
When the matter fields action is relativistic, we may write H ab in terms of just a scalar function of k, 

iU(fc) = n(«)K a <Ji(A), (20) 

where <5^,(fc) is the transverse tensor: <5^,(fc) = S ab — k a k b / n 2 . With this information, we can then produce 
more convenient expressions for S e ff(A), 

S eff (A) ~ i e 2 J -0^A:(k,O)n ab (k)A a (k,O) 
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7l F: b (k,0)U( K )F ab (k,0), (21) 



which shows that k 2 II plays a similar role to A for the real scalar field. 

When matter is nonrelativistic, the corresponding effective action, S^Jf, will have the same structure as 
in the first line of the previous equation, 

S?„(A) = \e 2 J -0- d Al{k,Q)Ull{k)A a {k,Q) (22) 



with an, in principle, different vacuum polarization tensor II™£. This tensor will still be transverse (current 
conservation) but it will fail to be just a scalar times the transverse 8. Lack of relativistic invariance means 
that its structure will be more complex, depending on two independent functions, which in turn will depend 
on |feo| and |k| separately (note that |k| is the norm of a d — 1-dimensional vector). 

Assuming rotation invariance on the (d — l)-dimcnsional spatial plane of each mirror, we can decompose 

n™£, using now two independent tensors: P^ b and P ab , that can be constructed having in mind that only 
the Galilean group of symmetries is relevant here. A convenient choice for those tensors is such that P^ is 
transverse in the d — 1-dimcnsional sense, namely, 

"00 = = -mm = "io ' P%j = "ij rrr j (23) 

where i, j denote run over the values 1, . . . , d — 1, and 

PL = Sab - Pab ■ (24) 

Note that, with the conventions above, they verify: 

(P x ) 2 = P- 1 , (PH) 2 = pH , P-LpH = , P" 1 - + PH = S ± , (25) 

and, on the other hand, each one of them is transverse to k a . 
Then the general structure of II™£ is: 

fis = n x p a \ + fin pI , (26) 

where TI 1 - and n" are functions of \ko\ and |k|. These two functions are related to different properties of the 
medium. In the T > case, for example, the function IP is responsible for the Debye screening of static 
electric fields, with a screening length determined by IP (ko — 0, |k| — » 0). Static magnetic field, on the other 
hand, are not screened, although there is non-vanishing screening when time-dependent magnetic fields are 
considered. They are also related to plasma oscillations, with plasma frequencies determined also by the 
coefficients above. 

Finally, to find the Casimir energy corresponding to the gauge field case, we recall that that energy can 
be derived from the Euclidean functional: 

Z = f [DA] e~ s{A ^ (27) 

where [DA\ denotes the integration measure including gauge fixing and, as in the scalar case, S is the sum 
of the free action plus the effective actions concentrated on the two mirrors: 

S eff = ^ 2 j dX * j dX 'l j ^^(Md)fiW« J 

2 

x S^ b (k)^2S(x d -a a )S(a a - x' d ) A b (k,x' d ) , (28) 

where we assumed relativistic matter fields. 

The path integral must be gauge-fixed; we adopt the 'axial' gauge A d = 0. This is a convenient choice, 
since it yields a quite simple form for the action: 

5 = 2 / dXd J dX ' d j (2^F Aa( - k ' Xd) Mab(k] Xd: X ' d) Ab ^' X ' d) ' (29) 



where 



M ab {k;x d ,x' d ) = -d 2 d 6(x d -x' d )6 l ab {k) 



Z 

+ { [(-dj + K 2 ) 6(x d - x' d ) + e 2 k 2 fi(«) [ £ 5(x d - a a )S(a a - x' d )] }# 6 (fc) (30) 
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and 5 ab = 8%, - 5 ab . 

The formal result for Z is 



Z 



_1 

(det[M ab (k;x d ,x' d )]) 2 , (31) 



where the determinant is evaluated over all the (continuum and discrete) indices. Since each variable in the 

path integral may be uniquely decomposed into transverse and longitudinal components: A a = A^ + A d , we 
have that: 

VA = VA 1 ^ VA} , (32) 

and, as a consequence, 

d-l _ 1 

Z = (detfM^A;;^,^)]) 2 (dct [MHk;x d ,x' d )}) * , (33) 

where 

M^{k;x d ,x' d ) = (-0 2 + k, 2 ) d(x d - x' d ) 

2 

+ e 2 k 2 U(k) [ 2J <5(^d - a Q )<5(a Q - x' d )] 

MllCfcjSd.si) = (-5 2 )<5(.T d -.^). (34) 

Note the power of d — 1 in (f33|) , which appears as a because 8^ is the identity on the space of transverse 
fields, ad-1 dimensional space. 

From the previous expressions, we see that Z is equivalent to the product of the vacuum functional for a 
massless free field which sees no mirrors, times d — 1 identical factors corresponding to scalar fields which 
interact with the walls with a (common) coupling X = e 2 k 2 II(k) . The free scalar field factor is irrelevant 
to the Casimir energy, since we subtract the energy when the mirrors are infinitely distant, while the others 
yields (d — 1) times a scalar field like energy 



Sof a ) = f ) - r dnn d - 1 ln(l - [l + 9 2 , J ~V 2 " a | 



(35) 



We specialize the previous study to an interesting concrete example: that of d = 3 with matter described 
by a massless Dirac field in a reducible An component representation (which preserves parity), we have |6j: 

(36) 

so that, after introducing the dimcnsionless variable x = kci , 



n («) = — , 

8k 




dablc a; = Ka , 
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$>(«) = g^s/ dxx 2 



Incidentally, this is one of the cases studied in [2|- This expression has the remarkable feature that the 
dependence of the Casimir energy with the distance is exactly equal to the one for a pair of perfect mirrors, 
albeit with a different global factor C{e 2 n) that interpolates between zero, for vanishing coupling constant 
e 2 = 0, and the result for perfect mirrors for e 2 — > oo. 

This is to be contrasted to the non-relativistic case which, using the same gauge fixing as before, yields a 
result composed of two contributions, corresponding to the modes that are spatially transverse or longitudi- 
nal: 



S ( a ) = „ * 2 il - r dKK d - 1 ln(l - [1 + J r 2 e' 2Ka \ 

0V ' 2<W 2 r(d/2) J X L e 2 K n^(«) J J 

/ dKK^ 1 \ a {l-[l + - T ^—]-\-^}. (38) 

Jo L e^Klfll(K ' 



2 d 7r d/2 TW 2) J Q I L e 2 K nll(K) 

Finally, we attempt to reinterpret the previous results for the Casimir energy as emerging from the existence 
of a discrete set of modes, reflecting the existence of (imperfect) boundary conditions at Xd = 0, a. We shall 
focus, for the sake of simplicity, on the case of a massless scalar vacuum field in d + 1 dimensions. After 
integrating out the matter degrees of freedom, one can evaluate the in-in, Schwinger-Keldysh, or closed time 
path effective action S^f F p}- On general grounds, this effective action has non-local kernels, denoting the 
existence of dissipative and noise effects. Taking the functional variation of the effective action, one can 
write a Langevin equation for the system field as 

U V + 9 2 Y. <tt'd d - l x\ G Tet (t,x h t\x\ l )S(x d -x a Mt',x\ l ,x a ) = J2 K ( t ^\\) S ( x d-^a)- (39) 

a a 

where G re t is the dissipation kernel, and K is a stochastic force whose correlation function is related to G re t 
via the fluctuation-dissipation relation. 

If we neglect the effects of the noise, the dissipative system will have complex eigenfrequencies. In order 
to compute them, we put the thin mirrors inside a larger box of length L, and impose periodic boundary 
conditions. Setting K = in Eq. (|39|) . it is possible to read the discontinuity of the field ip around the defect 
at Xd = and Xd = a as 

P'L+ -^L. =9 2 Y, [dt'd^x'u G rrt (t,X||,t / ,x||Mt',x(, > s a ) > (40) 

with a — 0, a. 

The condition that defines the eigenfrequencies is f(u>, a, L) = 0, with 

f(u, a,L) = \ (e 4 ^ L - e 2ik * a ) ( e - 2lk « a - l) - ^ (e uk « L - l) + ( p-) ~ (e 2ik * L - if = 0, (41) 

4 Arot V Aret / 



where kd = * oj 2 — k 2 and A rct = A re t(w, kn) is the Fourier transform of the retarded Green function G re t- 

Note that A rct is the appropriate analytic continuation of the Euclidean function A(k) to Minkowski space. 
Let us first consider the case in which A rct is real. When A re t is constant, it corresponds to a non dissipative 
situation, i.e. a free field with massive terms concentrated on the position of the mirrors [8|. On the other 
hand, when A re t depends on frequency, it is unrealistic to assume that it does not have an imaginary part, 
because of causality. In spite of this, one can compute the Casimir energy following a standard procedure. 
From Eq. (|41[) . it is easy to prove that, in the limit L —> oo, the admitted eigenfrequencies in the cavity are 
real. In this particular case the Casimir energy can be obtained as a sum-over-modes, 

£ = ^^2(^0,1 -Uoo,i), (42) 



which may be computed using the so-called generalized argument theorem 



^ dz Z — - 111 -77 ' - rr , (43) 

47^7 {2ttY~ 1 J dz f(z,L/2,L) y ' 

The integral in the complex plane reduces to an integral on the imaginary axis z — iy, and therefore, when 
A rct is an even function of u>, one gets 



«)4/ar* t .-.-^. + s] i. <«» 



2tt./ (2 7 r) d - i y " I Ar«t(*y,k|| 



which is consistent with the result obtained in Eq.JMJ) from the Euclidean effective action. 

However, when the effective strength of the coupling between field and mirrors comes from virtual effects 
of a matter field, one expects to have a complex spectrum. For example, for the case of massless Dirac fields 
in d = 3, the Euclidean function A is given by X(k) = ne 2 «/8 (see Eq. (|36|) ). Therefore, the retarded analytic 
continuation to Minkowski spacetime is |9| 



A ret (u;,k||) = -^- (i sign^)^ 2 - kjj0(o; 2 - kf ) - ^kf - w 2 0(kf - uj 2 )) , (45) 

introducing a dissipative term into Eq. (|39p . It is possible to show that, in this particular case, Eq. (|4"Tj) 
only admits propagating (o> 2 > k 2 ) complex solutions. In the limit L —> oo, they are given by 

9 /ran i , ,_ 16 A , , ,„. 

o& = In 1-— I +kf. (46) 

\ a a ne z J " 

where m is a natural number. 

As expected, in order to evaluate the Casimir interaction energy it is not correct to use the sum-over- 
modes formula Eq. (|42|) . since the modes of the effective theory have complex frequencies [10(. The correct 
procedure is to consider the sum over the real cigcnfrcqucncics of the full system, as implicitly done in the 
Euclidean approach. Alternatively, one should include the effects of noise in the effective theory, as done by 
Lifshitz in the derivation of the well known formula for the Casimir interaction in the presence of real media 
[111 ]. However, it is remarkable that the correct result derived using the Euclidean approach Eq.(fT4"|) formally 
coincides with the one derived from the sum over eigenfrequencies assuming an unphysical A rc t which is 
real and non constant Eq. (|44[) . A similar situation holds for Lifshitz formula, that has been derived using a 
sum over modes for unphysical permittivities with no imagina ry p arts [l2|, and also for the Van der Walls 
interaction between atoms immersed in an absorptive medium [13j. 

We conclude with some comments about the previous models, the assumptions made to motivate them, and 
the results obtained therefrom. An important point has to do with the relationship between the symmetries 
of the model and the resulting (approximate) boundary conditions. We have seen that, assuming rclativistic 
invariance and parity conservation on the mirrors' world volume, gauge invariance implies that the physics 
depends on only one function. This function, the scalar part of the vacuum polarization function, is only 
effective to impose a particular kind of boundary condition, which affects the components of the electric 
field that are parallel to the mirrors. Had one wanted to impose more general boundary conditions, one 
should have had to relax some assumptions regarding symmetries. This has been, indeed, explicitly shown 
for the case of a nonrelativistic system, where there appear two different functions. This situation also arises, 
for example, when the mirrors have a finite width. Finally, parity symmetry could have been broken, for 
example, by considering mirrors which are permeated by a strong magnetic field pointing along the direction 
of the Xd coordinate. This set up would produce exotic boundary conditions, mixing the transverse electric 
field with the normal magnetic field. 
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